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Geometric Brownian motion (GBM) is a model for systems as varied as financial instruments and 
populations. The statistical properties of GBM are complicated by non-ergodicity, which can lead 
to ensemble averages exhibiting exponential growth while time averages suffer collapse. A common 
tactic for dealing with this difference is diversification. In this letter we show that diversification 
will not solve the problem but only delay the collapse. 



PACS numbers: 02.50.Ey,05.10.Gg,05.20.Gg,05.40. Jc 

Geometric Brownian motion (GBM) is a useful model 
for systems in which the temporal evolution is strongly 
affected by relative fluctuations, such as stock portfolios 
and populations. Since fluctuations have a net-negative 
effect on growth in such systems, due to the multiplica- 
tive nature of the noise, diversification is a common strat- 
egy for assuring growth. In the language of statistical 
physics diversification involves a partial ensemble aver- 
age (PEA) over a finite number, N, of elements or tra- 
jectories generated by GBM. This raises important but, 
as yet, unanswered questions such as: How do the PEAs 
compare to the ensemble average (N — > oo)? If there are 
significant differences how and when do they arise? In 
this letter we analyze PEA's of the GBM both analyti- 
cally and numerically. Our main result is that diversifica- 
tion is a temporary solution. In our examples diversified 
"portfolios" only exhibit growth for a short time before 
arbitrarily small fluctuations lead to collapse. In other 
words the diversified portfolios are linearly unstable. 

Geometric Brownian motion is defined by 

dx = x(fidt + adW), (1) 

where /i is a drift term, a is a noise amplitude, and 
W(t) — J Q dW is a Wiener process. Without the noise, 
i.e. a — 0, the model is simply exponential growth at 
rate fx. With a ^ it can be interpreted as exponential 
growth with a fluctuating growth rate. 

To solve (Eq. [I]), one computes the increment dln(x) 
of the logarithm of x, integrates and exponentiates. The 
interesting step is computing the increment because this 
requires stochastic calculus. In writing (Eq. [I]) we had 
in mind an interpretation of the equation in the Ito con- 
vention pp. With this convention, it is well known that 
dln(x) = (p — ^-)dt + adW, which by exponentiation 
implies the solution 

x(t)=x(0)expUfi-^)t + aW(t)j. (2) 

For simplicity, we will assume the initial condition x(0) — 
1. Figure [T] illustrates the nature of this process. 
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FIG. 1: Percentiles from 95 to 5 (top to bottom) in steps of 
5, based on 10,000 realizations of x(t). The parameters (used 
for all illustrations here) are = 0.05 and a 2 = 0.2. The red 
straight lines show the ensemble average (upward sloping) and 
an exponential decreasing with the time-average growth rate. 
The ensemble average is essentially meaningless for a single 
realization, whereas the time-average growth rate accurately 
describes the typical behavior. 

The process x(t) is not stationary. This implies that 
where averages can be defined, there is no guarantee for 
ergodicity, i.e. the equality of ensemble and time aver- 
ages P]. The time average of the process itself is (a.s.) 
either or diverges positively, whereas the ensemble av- 
erage is an exponential function of time. To capture the 
non-ergodicity of the process in well-defined averages of 
an observable, we define the following estimator for the 
exponential growth rate 

Sert(t,JV):=^In«a:i(t)M, (3) 

where we call (■) N := jf'^2 i • a PEA. The estimator 
looks at the growth rate of a PEA (it is not a PEA of the 
growth rate) , i. e. the logarithm is taken outside the aver- 
age. This is crucial to leave the non-ergodic properties of 
the process, x(t), intact. The time-average growth rate, 
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denoted g, is found by letting time remove the stochas- 
ticity in the process. Mathematically, this is the limit 



g := lim g cst (t,N) 
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(4) 



The ensemble-average growth rate, denoted (g) , is found 
by letting an increasing ensemble size remove the stochas- 
ticity. Mathematically, this is the limit 



(g) := lim g cst (t,N) = fj,. 
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(5) 



The non-ergodicity of the process is manifested in the 
non-commutativity of the limits lim t _ i . 00 and lim7v_ i . 00 . 

Both ensemble and time averages are mathematical ob- 
jects and therefore separated from physical reality by the 
divide that separates logic from matter. Nonetheless, 
both averages carry practically meaningful messages. To 
identify the regimes where they "apply", that is, where 
they reflect typical behavior, it is important to under- 
stand more about the general case where both the obser- 
vation time, i, and the ensemble size, N, are hnite and 
arbitrary. 

In [3J, the time-average growth rate, (Eq. |4|, was com- 
puted for a single system, N = 1, by letting t diverge. 
But the case of arbitrary N was not treated. The ensem- 
ble average was computed for arbitrary t. But the limit 
N — > oo was not taken explicitly, relying on the fact that 
in this limit the PEA, (xi(t)) N , is the expectation value, 
(x(t)}. Below we show that (Eq. [4| holds for arbitrary N 
and characterize the process of the convergence of (Eq. [5| 
as N — > oo. 

We begin by showing that for a single instance, N = 
1, the distribution of g cs t(t, N — 1) approaches a delta 
function centered on fi — a 1 /2 in the limit t — > oo. 

Substituting (Eq. |} in (Eq. [3}, g est (t,N=l) = fi - 
^ + \aW(t). We know that the distribution of W(t) is 
Gaussian with mean and standard deviation i 1 / 2 , which 
we write as P(W(t)) = Af(W(t); 0, t 1 ' 2 ). To compute 
the distribution of g cs t(t, N=l), we use the transforma- 
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tion law of probabilities, P(g) — P(W) 
dg/dW = f and solving (Eq.^j) for W(g), this yields 




P(g es t( t > N = 1 ))=N Sest;M' 



The limiting behavior of this distribution for t 
the Dirac delta function 



lim P(g cst (t,N=l)) = S g cst - ( M 
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In other words as t — > oo, the observed growth rate will 

2 

differ from fj, — with probability zero. 

Next, we consider N instances of (Eq. [l]). At each mo- 
ment in time, the TV instances are averaged, as illustrated 



(x(t)>. 
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FIG. 2: Typical trajectories {x(t)) N for different ensemble 
sizes N, from top to bottom TV = 1024, 512... 1. For each 
ensemble size, 10,000 trajectories were generated. At each 
time the median, the 49th and 51st percentiles are shown. 
The solid straight red lines illustrate the ensemble-average 
and time-average growth rates. The single system (bottom) 
is dominated by time-average behavior; for increasing TV the 
trajectories stay close to the ensemble average. 



m Fig.|2j and we are interested in the long-time behavior, 
t — > oo, of the object (x(t))»j- 

Again, substituting (Eq. in (Eq.|3|, 

exp ( (fi - °^]t + aWi(t) 
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g est (t,N) = - In 



N 
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In ((exp (aW t (t))) N ) 



(8) 
(9) 
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The difficulty with this equation is the logarithm of an 
average of exponentials. Unlike in the case of the single 
system, the logarithm does not simply undo the exponen- 
tial, and the non-trivial behavior of typical trajectories 
of PEAs is a direct result. 

To show that (Eq. [4]) holds for arbitrary TV, we will 
proceed in two steps by showing that as t — > oo 



• the probabilty of finding g est (t,N) > /i 
proaches zero (upper bound) 



ap- 



• the probabilty of finding g est (t,N) < \i — ~ ap- 
proaches zero (lower bound). 

Upper bound: 

Equation (|9| is a growth rate estimate of an average. But 
the average cannot be larger than the largest individual 
term. This establishes an inequality, namely an upper 
bound on g est (t,N). 



y 



-a max Wi(t). 
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g cs t(t,N) < fx 

A value g es t(t, N) > /i — ^- + e is thus only possible if 
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maxWj(£) > — . 
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(11) 
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The probability of such an extremum is [3] 
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WO; 0, Vt)dz 
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Two interesting properties can be observed. First, for 
finite N, and e > 0, 

lim P f max Wi(t) > - ) = 0, (13) 
t-*oo \ * cr I 

the desired result. This is because the width of the dis- 
tribution J\f(x\ 0, y/i) increases as y/i, whereas the upper 
limit of the integral increases as t, outpacing the diver- 
gence of the width. In the limit t —¥ oo, the entire distri- 
bution is integrated, which yields 1 due to normalization. 
Second, for finite t, the limit 
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lim P ^niax W{ (t) > 
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This must be so because the term that is being raised 
to the 7V th power is less than 1 for finite t and therefore 
vanishes exponentially with N . In other words, in the 
limit of diverging ensemble size, the method fails to give 
an upper bound on the estimated average growth rate. 
This is so because maxf Wi{t) diverges in this limit. 
Lower bound: 

The lower bound on g cs t(t,N) is obtained in the same 
way as the upper bound, by switching the inequality and 
considering the minimum of the N instances at time t. 

We have shown that as t —> oo the probability for 
observing a growth rate g es t{t,N) ^ fj, — ~ approaches 
zero for any finite N. 

In proving (Eq. El), we have used extreme values. This 
is no coincidence: extreme values will be the key to un- 
derstanding our problem. This is because the exponential 
introduces a weighting that leads to a finite contribution 
to the average from extreme values whose relative fre- 
quencies vanish in the limit N — > oo. Considering the 
discrete-time, discrete-space random walk, it is clear that 
the absolute maximum - not the typical maximum but 
the largest possible value - scales in a light-cone fashion 
as t and not as t 1 / 2 . This is reflected in two regimes of 
actual physical diffusion, a short-time regime where ex- 
trema among the positions of diffusing particles scale as 
t and a long-time regime where they scale as t 1 / 2 , beau- 
tifully illustrated in [5] . Further illustration of this argu- 
ment using the discrete multiplicative binomial process 
can be found in [BJ. For the Wiener process, we cannot 
consider the largest possible values because they are infi- 
nite for any t > 0. This is a well-known limitation of the 
model. As pointed out in [7], the canonical solution to 
the diffusion equation violates special relativity because 
it allows diffusing matter to exceed the speed of light. 
This is visible in the small-i behavior: while for large 



t positions at a distance y/i correspond to slow motion 
and the Wiener process is a good model for the physics 
here, at short times the speed of such motion diverges 
and becomes unphysical. 

In [BJ it was argued that in order to observe ensemble- 
average behavior (N — > oo) for a time r in a PEA, 
N ~ cxp(r) multiplicative systems are required. This 
scaling is a consequence of the fact that the probability 
of t consecutive up-moves in a random walk decreases 
exponentially with r. The multiplicative nature of the 
process enhances large outliers and leads to the extreme 
values dominating the (linear) average behavior. After a 
time r ~ ln(iV) the absolute extremes become a-typical 
for the ensemble size, and their effect diminishes, lead- 
ing to a deviation from ensemble-average behavior. It is 
important to note that the result in [BJ is derived in the 
large N limit. Here we reconsider this problem, phrasing 
it in terms of the stability of PEAs, and obtain a some- 
what different conclusion. In particular we find that the 
PEA deviates from the ensemble average at an earlier 
time and is linearly unstable, i.e. unstable with respect 
to arbitrarily small perturbations coming from the noise. 

We begin by defining the deviation ejv(f) of the PEA 
from the ensemble average by 



(x(t)) N = exp(iit) + e N (t). 



(15) 



Initially, trajectories will approximate those of the en- 
semble average so that we can approximate the deviation 

as 

e N (t) = exp^fd + a^ dW^j ^ - exp(/it). (16) 



Replacing 



Jo" dw ) N b y \J((Wi(t)) 



jf we ob- 



tain an expression for the scaling behavior of the devia- 
tion 



ejv(i) ~ cr exp(^i) 



(17) 



It can be shown that ejv(i) is the PEA of the solution to 
d,CN=i(t) — dtfi£N = i(t) + acxp(nt)dW (18) 
with the initial condition e^r(i) = at t = 0. In addi- 



tion equation (Eq. 17 1 is the lowest-order contribution 



from an asymptotic series generated from an 
iterative solution to (Eq. [lj (manuscript in preparation) . 

The approximation in (Eq. |16[ ) neglects any non-linear 
effects. Nonetheless it is informative of the trade-off be- 



tween N and t. We can set (Eq. 17) equal to some finite 



value and derive an expression for the time, r, it takes to 
reach this deviation for a given ensemble size N . From 
(Eq.[l7|, this is 



- (ln(ejv(f) 



In (a a 



(19) 
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Compared to the logarithmic scaling, r ~ ln(iV), (Eq. 19 ) 



includes a correction (the second term on the RHS). For 
large characteristic times, r 3> | In (cr-y/r) \ (i.e. large 
values of eN{t)y/N), we can neglect this correction. Note 
however, that for the asymptotic expansion to be a valid 
approximation we must have a^/t/N is small for t = r. 

In Fig. [3] we show the times r where absolute deviations 
of magnitude 0.1 and 1 of the PEA from the ensemble 
average were reached in the trajectories in Fig. [2j 

Considering the asymptotic series from the iterative 
solution to (Eq. [I]), we note that the ratio of the first- 
order term to the second-order correction for the data 
in Fig. [3] is less than 0.55 for N > 10. For e N (t) = 1 
where the ln(iV)-scaling appears to be valid, this ratio is 
small, namely ranging from 0.55 for N = 10 to 0.08 for 
N = 1000. The ratio of the contribution of the second 
to the first order term is even smaller when ejv'(t) =0.1. 
This implies that the linear approximation is a good de- 
scription for a wide range of parameters. 
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FIG. 3: Time where the median for a given N, shown in 
Fig. [5] first differs from the ensemble average exp(fit) by at 
least 1.0 (ti, top) and 0.1 (to.i, middle). Error bars show 
the corresponding times for the 49th and 51st percentile. A 
purely logarithmic fit of the form r = a+b ln(iV) works well for 
deviation 1.0 (with a — 9.46 and b = 3.88). For deviation 0.1 
occurring at smaller t corrections become visible; to compare 
shapes the solution for t(N) of (Eq. 19 1 is shown in arbitrary 
units. Bottom: the percentage difference, A, between the top 
data and their fit. 



Also features associated with much larger deviations, 
such as zero-crossings of the growth rate, (Eq.[3]), approx- 
imately follow logarithmic scaling. This can be seen in 
Fig. |2j where the spacing of successive zero-crossings of 
the growth rate (i.e. trajectories crossing 1) is approx- 
imately constant for each doubling of N. For a small 
deviation ejv(i) = 0.1 we find a non- logarithmic shape 
similar to (Eq. |l9|) including the correction. 



The fact that properties of a linear approximation feed 
through to the full non-linear solution is surprising but 
not unheard-of. 



Equation (Eq. 18 1 is identical to the Cahn-Hilliard- 



Cook theory which describes the evolution of a class of 
materials after a quench into an unstable state 8J. The 
implication is that growth associated with systems mod- 
eled by PEAs of GBMs such as populations, stock port- 
folios and the stock market itself is inherently unstable. 
This leads to the conclusion that any PEA will eventually 
be dominated by the same time-average behavior ((Eq.Q 
holds for arbitrary N) . These results have important im- 
plications for the relevance of diversification strategies 
under realistic conditions in fields ranging from financial 
risk management to ecology and evolutionary biology. 
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